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Atomic wave interferometers are tied to a telescope pointing towards a faraway star in a nearly 
free falling satellite. Such a device is sensitive to the acceleration and the rotation relatively to the 
local inertial frame and to the tidal gravitational effects too. 

We calculate the rotation of the telescope due to the aberration and the deflection of the light 
in the gravitational field of a central mass (the Earth and Jupiter). Within the framework of a 
general parametrized description of the problem, we discuss the contributions which must be taken 
into account in order to observe the Lense-Thirring effect. 

■^Jsing a geometrical model, we consider some perturbations to the idealized device and we calculate 
tH^corresponding effect on the periodic components of the signal. 

GJbme improvements in the knowledge of the gravitational field are still necessary as well as an 
idciilease of the experimental capabilities ; however our conclusions support a reasonable optimism 
fos-the future. 

Anally we put forward the necessity of a more complete, realistic and powerful model in order to 
obtain a definitive conclusion on the feasibility of the experiment as far as the observation of the 
JfSRpe-Thirring effect is involved. 
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00 ' Clocks, accelerometers and gyroscopes based on cold atom interferome- 
. try are already among the best which have been constructed until now and 
further improvements are still expected. With the increase of the experi- 
mental capabilities it becomes necessary to consider more and more small 
effects in order to account for the signal, therefore (relativistic) gravitation 
has to be considered in any highly sensitive experiments, no matter what 
O they are designed for. 

The performances of laser cooled atomic devices is limited on Earth by 
5h 1 gravity. Further improvements demand now that new experiments take 
^P" place in free falling (or nearly free falling) satellites. A laser cooled atomic 
J> , clock, named PHARAO, will be a part of ACES (Atomic Clock Ensemble in 
Space), an ESA mission on the ISS. Various other experimental possibilities 
involving " Hyper- precision cold atom interferometry in space" are presently 
considered. They might result in a project (called "Hyper") in the future 
' 

The aim of the present paper is to hold the bookkeeping of the various 
gravito-inertial effects in a nearly free falling satellite. For this purpose we 
consider the most ambitious goal which has been considered for Hyper i. e. 
the measurement of the Lense-Thirring effect. 
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The Lense-Thirring effect is a local rotation of a gyroscope relatively to a 
telescope pointing towards a far away star. It is a relativistic consequence 
of the diurnal rotation of the Earth which "drags the inertial frames" in 
its neighborhood. 

The angular velocity of the telescope relative to the gyroscopes depends 
on the position. Therefore, in a satellite, it is a function of the time. In 
Hyper, the angular velocity is measured by atomic-wave-gyroscopes and 
its time dependence is analyzed |9j- The consequence is that the device is 
sensitive to the variation of the gravitation in the satellite and no to the 
gravitation itself. We do not believe that it is easy to achieve the required 
stabilization of the gravitational field due to the local masses but it is not 
impossible in principle. For this reason we will study only the tidal field of 
far away masses whose effect cannot be removed at all. 

The parameter which plays a role in the calculation of the Lense-Thirring 
effect is the angular momentum of the central mass. It is much bigger for 
Jupiter than the Earth. Therefore we will discuss both cases, without any 
consideration on the cost of the corresponding missions. 

In the sequel the greek indices run from to 3 and the Latin indices 
from 1 to 3. We use the summation rule of repeated indices (one up and 
one down). 

The Minkowski tensor is r] a p = diag [1, —1, —1, —1]; its inverse is rj a ^ . 
We use geometrical units where c = G = 1. 

II. THE LOCAL EXPERIMENT IN A SATELLITE 

In the satellite, the experimental set-up consists in a telescope pointing 
towards a far away star in the 'u'(i) direction and two orthogonal atomic 
Sagnac units in the planes [^(3) ; (i)] an d ["^(2)>^(i)l °f n g-ffl 
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FIG. 1: The experimental setup. 
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A. The atomic Sagnac unit 



An atomic Sagnac unit (ASU) is made of two counter-propagating atom 
interferometers which discriminate between rotation and acceleration (see 
figure |2[a). 
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FIG. 2: An atomic Sagnac unit (ASU). 



Each interferometer is a so-called Ramsey-Borde interferometer with a 
Mach-Zehnder geometry ( figure 0-6). The atomic beam from a magneto- 
optical trap interacts three times with a laser field. In the first interaction 
zone the atomic beam is split coherently, by a Raman effect, into two beams 
which are redirected and recombined in the second and the third interaction 
zone. 

The mass of the atom depends on its internal state, therefore it is not 
a constant along the different paths. However, the change of the mass is 
very small; it leads to negligible corrections on the main effects which is 
already very small. In the case of the cesium, the mass is m = 133 x 1.66 x 
10~ 27 = 2.2 x 10~ 25 kg and the wave length of the lasers is A = 850nm. 

Airh 

The momentum transferred to the atom during the interaction is — — . The 

A 

recoil of the atom results in a Sagnac loop which permits to measure the 
angular velocity of the set-up relatively to a local inertial frame. The device 
is also sensitive to the accelerations. 

In an ideal set-up the two interferometers are identical coplanar paral- 
lelograms with their center Os and 0' s at the same point but many per- 
turbations have to be considered. The geometry of the device is actually 
determined by the interaction between the initial atomic beam and the 
lasers ; Therefore a full treatment of the atom-laser interaction in a grav- 
itational field is obviously necessary to study the response of the Atomic 
Sagnac Unit (ASU). However the geometrical model is useful to give a 
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physical intuition of the phenomena. In this context we assume that the 
two interferometers remain idealized identical parallelograms but that Og 
and O'g are no longer at the same point : This is the only perturbation 
that we consider here. It is sufficient to take the flavor of the gravitational 
perturbations which have to be taken into account and, more generally, of 
the difficulty inherent to such an experiment. 

B. The phase difference 

Let us assume that the fundamental element is known in some coordi- 
nates comoving with the experimental set-up : 

ds 2 = (l + K (om )dT 2 + 2K mk) dTcLXW (1) 
+ {V( km +K (km )dXWdXM 

In order to calculate up to first order the gravitational perturbation of 
the phase due to Ki a )(p), we use a method which we summarize now 0- 
First we calculate the quantity \& : 

* = K {0){0) + 2K mk) v^ + K [m) v^v<j) (2) 

(k) 

where v g is the velocity of the atoms [i.e. the unperturbed group veloc- 

ity). 

The quantity ^ is a function of the time and the position of the atom. 

Then we consider an atom which arrives at time t at point B of figure Gi- 
ft . Now the position is a function of the time t' only because t is considered 
as a given quantity. The function 4" is a function of the time, t' , only. The 
phase difference is 

* (0 dt '~^ * (*') dt ' ( 3 ) 

-2T D /(2) 2 Jt-2T D /{1) 

The integrals are performed along path (2) and (1) of the interferometer 

9 

7TL C 

(figure The "angular frequency" ui is defined as — — . 



C. The local metric 

In order to calculate 8(p, we must know the local metric Gt a \i^\ = 
V{ a )(i3) +K [a){fj) . 

We choose an origin, O, in the satellite and at point O a tetrad 

| u (o)5^(j.)) u ^))Wpjl where is the 4- velocity of point O and where 
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the three vectors = {^^{k)} are represented on the figure^ The 

vectors of the tetrad are orthogonal : UuJi a ( ff ) = f]{p){cr)- The coordinate 
indices, a, /3, a, etc, are lowered or raised by the means of the metric tensor, 
<7a/3 or g a @ ; the Minkowski indices are raised or lowered with the Minkowski 
metric or its inverse, rj a p or rj a ^ . An Einstein indice, a, can be changed 
into a Minkowski indice (p), by the means of the tetrad and vice versa : 

«fo()a = ()(,) and «(p)() (P) = Q - 

The tetrad is the natural basis at point O of comoving coordinates X^ a \ 
The proper time at the origin is X^°\ The space coordinates are the X {k \ 

Once the origin and the tetrad are chosen, the metric at jDoint M and 
time t is expanded relatively to the space coordinates of M 



ds 2 = G [a)([j) dX {a Ux {l3) with (4) 



G (0)(0) 



l + 2a-I+(a.i) 2 -(sixl) 2 - #(o)(fc)(o)(j) 
- l -R mk)mjU t)X {k) X^X^+... (5) 

G(0)(m) = 0(m)(fe) - g-R(0)(fc)(m)(j) X^ k > 

-\R mmm)Uue) X^X^X^ + .. 

G( n )(m) = V(n)(m) - g -R(n) (fc) (m) (j) X ^ X ^ 

-QR(n){k)( m )(j)M) Xik)Xb)xW + ■- 

where we have used vector notations i.e. ~a for {a^- 1 }, ~a ■ X for 
^o"lW, etc. Every quantity, except the space coordinates are 
calculated at point O. Thus they are functions of the time T = X^°\ 
R(a)(/3)(<r)(fj,) is the Riemann tensor obtained from R a f3 a ^ at point O : 

Retftafi — ^a—f3fi : a ^a—j3(T,/j ~t~ ^ e — afi fin^ e—cta (^) 

where r a _^ iff is the Christoffel symbol. 
^(j)(fc) is the antisymmetric quantity 

Q (Mk) = 2 (-9(0)0), (fc) -5(0)(fe),(j)) o (7) 

^ x x 




o 



6 



Due to the antisymmetry of fl( m )(k), the quantity fi( m )(fc) X^dX^ 
which is present in the expression of ds 2 can be written as 
n (m)(k) X^dX^ = (Jl A ■ dX. The space vector Tt is the physi- 
cal angular velocity. It is measured by gyroscopes tied to the three space 
orthonormal vectors u?m : 

The vector ~a is the physical acceleration which can be measured by an 
accelerometer comoving with O. It is the spatial projection at point O of 
the 4-acceleration of point O. 

At point O {i.e. X = ) the time T is the proper time delivered by an 
ideal clock comoving with O. 



III. FROM THE GEOCENTRIC COORDINATES TO THE 
COMOVING COORDINATES 



A. The geocentric coordinates 

We define the time coordinate x° — ct and the space coordinates x k . 
We use the notations r — {x k } — {x,y,z} and we define the spherical 
coordinates {r, 9, ip} , i.e. x = r sin 9 cos ip , y = r sin 9 sin ip , z = r cos 6. 

We consider a satellite and a point O which is the origin of the local 
coordinates in the satellite. We assume that the position of O is given 
by its three space coordinates, r — {x,y,z} = {x k } , considered as three 
known functions of the coordinate time, t. Then we define the velocity of 
, dr 

point O as v = — . 
1 dt 

The proper time at point O is s = T = X^K The motion of O can be 

described as well by the four functions x a = x a (s) . The four-velocity is 

dx a 

defined as u a = — — . 

ds 

In the sequel we consider the Parametrized Post Newtonian theories 
P|. The relevant PPN parameters which appear below are 7 and a%. The 
parameter 7 is the usual parameter connected to the deflection of a light 
ray by a central mass. The parameter a.\ couples the metric to the speed, 
—w, of the preferred frame (if any) relatively to the geocentric frame. In 
general relativity, a\ = and 7=1. 

Let us define now several quantities which will be used in the sequel : 

• 2M® is the Schwarzschild's radius of the central body (i.e. the 
Earth or Jupiter). As we use geometrical units, M@ is also its "mass". 

• J@ is the angular momentum of the central body in geometrical 

1 + 7 + a!i/4 -> 

units. The relevant quantity which appears below, is J = J®. 



We define J = 



J® — J® 
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J A r 1 

• g = — 2 — - 1 — otiU w is the definition of g, where w is the 

r J 2 

velocity of the observer, relative to the preferred frame (if any). 

• is the Newtonian potential 

^ = ^(W 2 (^) 2 P 2 + A) + 14 (8) 

where i?@ is the radius of the central body and [/* the potential due to its 
satellites, the Sun and the planets [Tof. In spherical coordinates the Legen- 

dre polynomial P 2 reads Pi — — (3 cos 2 9 — l) . The quadrupole coefficient 

is J 2 and A represents the higher harmonics. It depends on the angle ip 
and on the time t because of the rotation of the central body. 

In the non rotating geocentric coordinates the significant fundamental 
element is 

ds 2 = (1 - 2U) dt 2 + 2g Qk dx k dt - (1 + 2 7 C/) 5 jk dx j dx k (9) 

where {g) k — — (g) k — gok- In eq. ©, we have dropped post Newtonian 
corrections which are too small to be considered here. 



B. Orders of magnitude 

Table 1 below gives the order of magnitude of the various parameters 
which have been introduced previously. 





M® 


J© 


R® 


J2 


A 


Earth 


4.4mm 


145cm 2 


6400km 


- 10- 3 


- io- B 


Jupiter 


1.4m 


1700m 2 


71300km 


- 10~ 2 


< 10^ 



Table 1. 



In order to describe the physical situation we introduce four parame- 
ters : £, E, r\ and ji. 

First we define the order of magnitude 0\ — . / ® . The quantity (Oi)™ 

v R ® 

is denoted by O n . 

Then we consider a nearly free falling satellite on a nearly circular orbit 
of radius r ~ R®/£. This expression gives the definition of £. The velocity 
of the satellite is of order v = ^^ 2 Ox [Tl| . 

Now we define d = R@Oi and e such as X = ed where X is the size of 
the laboratory. 

We define rj. The velocity of the atoms is v g = r}0\. 
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Finally we assume that the various quantities such as the position of O or 
the geometry of the experimental set-up is known with a relative accuracy 
of order of /i. 





relative accuracy : [i 


Orbital parameters 


set-up parameters 


radius r = — ® 
velocity v = Oi^ 1 ' 2 

period T= f i/ 20l ce 


size X = R@Oi e ~ 60cm 
atom velocity v g — r\0\ ~ 20cms~ 1 /c 

Drift time 2T D = X/v g = R® ^ ~ 3s 



Table 2. : definition of Oi,£,£ and 77 



With £ ~ 0.9 one finds 





Ox 


e 


V 


r 


V 


T 


Earth 


2.6 10" b 


3.6 10~ 3 


2.7 10~ b 


7000km 


2.9 10~ b 


5900s 


Jupiter 


1.4 10~ 4 


6.0 10~ D 


4.8 10~ b 


78400km 


1.5 10"° 


12300s 



Table 3. 



C. Comoving non rotating coordinates 

We consider the following tetrad, e?, comoving with O : 



el = u ° = l + ^+U + 4 , e k =u k = (l + ^ + U) v k + 
e\ = (l + ^ + u) v k + 1 Uv k - g ok + £ 2 4 (10) 



4 - 5{+ l -^v k + l - 1 U8{ + eO i 

The local metric is derived from 0] with the change in the notations 
(a) — > d and it" CT j — > ej?. 

We limit the expansion of the metric at order e 2 £ 3 / 2 Og; therefore we 
consider only the linear expression of the Riemann tensor (eq. ©) and we 

c 2 

assume that the free fall is under control : ||<z|| << eO^I 2 x — (i.e. 

X 

\\a\\ « 8ms~ 2 for the Earth, ||a|| << 21 ms for Jupiter which is not very 
restrictive) therefore we neglect the term ■ X^j in the metric (JJJ. One 
finds 
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G M = l + ^a-X-U, yj X~ k Xi -^U, kfl X k XlX l + e 2 e0 6 (11) 

_ -\ rh 



{«o AXj + e 2 e /2 5 and G, Wl = %A + e 2 £ 2 4 (12) 



where fio is given below (see eq. Ulcty l while the expressions such as 

are nothing but (u, mn e ^ e j ) • The position of the observer changes with 

time, therefore this quantity is a function of T. 

We did not consider the time dependence of the potential U. One can 

AU r AU . 

prove that it is correct when x << £C?2 where is the relative 

U cT c U 

change of the potential during the time T c , at the distance r of the origin. 
This is generally the case. 

In G50, ^ ne accuracy is limited to the terms of order of e 2 ^Oq. One 
can check that in such a case, the approximation e™ = <5|™ is valid there- 
fore (U,kj) and U,fet= {U,kjt) - The same holds true for fig 

k / \ k 

fio) (see eq. J130 ). Therefore, one can identify the space 



(60) 



vectors e t of the tetrad and the space vectors d k of the natural basis 
associated to the geocentric coordinates. This would not be valid with an 
higher accuracy where terms smaller than e 2 £ 3 Og are considered. 
Calculating Qq one finds the usual following expression Q 



* ( J 3 



(j.n) 



(Q ds y ~ ((1+7)TOA«) 1 and (n Th )' 



(13) 
(14) 



1 -. dv\ 
r A *) (15) 



n is the direction of the satellite (fig. El, ^lt is the Lense-Thirring angular 
velocity, fids and &Th are the de Sitter and the Thomas terms [T^ : 



Table 4. 



With £ ~ 0.9, one finds 
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J®/ M i> 


n LT 




Earth 


750 


-lO^rads- 1 


^lO-^rads- 1 


Jupiter 


855 


-lO^rads- 1 


^lO-^ads- 1 



Table 5. 



D. Aberration and deflection of the light 



In the satellite, the experimental set-up is tied to a telescope which 
points towards a "hxed" star (see figure We assume that the star is 
far enough for the parallax to be negligible. However the light rays suffer 
a gravitational deflection from the central body and an aberration which 
depends on the position and the velocity of the satellite. These effects 
result in an angular apparent velocity which must be compared to the 
Lense-Thirring effect. 




FIG. 3: The deflection of the light. 



In space time, the direction of the light from the star is given by the 

{dk f ] 
1, — — > where ip is the phase of the light. In order to 
do <p J 

calculate the phase ip(t,x k ) at poin t \x k } and time t we use the method 
which is summarized in paragraph III Bl Now the line element is given 
by eq. JO}, and u> is the angular frequency of the light at infinity. These 
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calculations are developed in another publication Here, we just give 
the useful results. 

The main gravitational contribution is due to the monopolar term of the 
Newtonian potential : 

L a = \l, -^ + (l + 7 )— ^— —..+ 5t k + 5L k ) (16) 

where I is the unitary vector of figure |31 and ~n — ~r jr . 

• The term S£ k is due to the quadrupolar term of the central mass. 
This term is of order J2O2 when £ is nearly orthogonal to the plane of the 
orbit. 

• The term due to — ct\U w , a part of gok in the metric @, results 

(OL\ W ■ I \ 
1 . 
2(1+7)/ 

• The contribution due to the rotation of the central body is of order 
of J@/r 2 . The corresponding angular velocity is of order of J@/r 2 /T ~ 

— % — — ~ Qlt « Qlt- It is negligible. The same conclu- 

M| 2tt R® 2tt 

sion holds for the term A in eq. © ■ 

• The Sun, the satellites and the other planets, give a contribution 
due to £/» in (JSJ; it varies slowly with the time and it is negligible, especially 
within the framework of a Fourier analysis at a much higher frequency. 
An exception concerns the two satellites of Jupiter, Andrastea and Metis 
whose period is approximately 25 10 3 s which is the order of the period of a 
satellite on a low orbit. However their mass do not exceed 10 17 kg and the 
gravitational deflections remain completely negligible. 

For the observer O, the space direction of the light is the four vector 
A" = L a — Lpu 13 u a . The components of A" relatively to the tetrad are 



{a 5 H(o, a) 



We define A = A A such as -A a A Q = A ■ A = 1. 

The tetrad (|10fl is especially useful to catch the orders of magnitude of 
the various terms involved. However it is not the comoving tetrad that 
we are looking for because the telescope that points towards the far away 
star rotates relatively to this tetrad. The angular velocity of the telescope 
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relatively to |e|| is 0» = A A Straightforward calculations give 



3 (* 



dv 
v A — 
dt 



(17) 



M 1 + 7 
r 2 1 - n ■ I 



(lA vj k + (Ia nj 



l-v - 



l-n-t 



Let us notice that we neglect the terms of order - x £ 2 04, which are much 

r 



smaller than the Lense-Thirring angular velocity because — S- 

M ® 



>> 1 and 



E. Local coordinates tied to the telescope 



Now we introduce the tetrad tied to the telescope and the interferometer 
u9 a y It is obtained from e> through a pure space rotation (i.e. u? \ = e? = 
a") and whose vector u?U points towards the far away star = — A Q ). 

At the required accuracy, it is possible to give a description of the Hyper 
project with the Newtonian concept of space. 

The rotation of the tetrad relatively to | e ^| is characterized 

by the most general angular velocity fl u / e — fi* — wA where — zuA is an 
arbitrary angular velocity around the apparent direction of the star. The 
change of the tetrad u?^ < — ► ej? is just an ordinary change of basis in 
the space of the observer O. In this transformation, dT, Gqq = G(o)(o)> 
G^dX™ = G mk) dX^ and G Arn dX" l dX fl = G (j)(k) dX^X^ behave 
as scalars. We obtain the local metric from eqs. JTJ and (|12|) . Then, using 
the expression |21 of a straight forward calculation gives : 



* = 2d-X-U, [k){3) xWxW-±U, (km(e) X^X^X^ (18) 



!^{(l!fl + !J,)A^ v, 

(k) 



13 



n 



LT 



£A 



dv 
It 



AX 



V g + 



with |^O + ^*)ax| • ~v g = 

V£ e /2 o 5 . 

The Lense-Thirring contribution to * is ^ lt ~ --jj^etiOq. Therefore 
within the present framework, the expected accuracy is of order of 





Earth 


Jupiter 




18% 


1.5% 


( J©/M|) eevOe (J©/M|) V 



Table 6. 



IV. THE PHASE SHIFT 



Let us assume that any quantity can be known with an accuracy fi ~ 
1CP 4 . This condition is not restrictive for the orbital parameters and does 
not seem out of the present possibilities as far as the geometry of the 
experimental device. 

We consider that ^ is the amount of two terms, Vf^ and ^ u : the term 
is known; it can be modelled with the required accuracy while is 
unknown. The terms fulfills the condition /i x $t < £ 2 £?0§. With the 
previous orders of magnitude one finds 







Earth 


Jupiter 


e 2 eo 6 


eo 2 ~ 


1.6 10 5 e * u 


5.7 10 3 e tf u 






1.5 10~ 2 e * fe 


4.5 1CT 5 e * fe 


e 2 eo 6 






m 


1.2 10 3 e * u 


4.5 10 2 e tf u 




e$0 2 


e 2 eo 6 


m 

e£ 1 /2 0l n 


3 10- 2 e * fe 


6 10" 2 e * fe 



Table 7. 



ty u reads 



* u = 2 a ■ X - U, mj) I^I* 3 ' - 2 {fi A x} ■ V s + e 2 £ 3 6 /m (19) 

where the contribution U,(k)(j) X^X^ needs to be defined with an accu- 
racy better than e 2 £ 3 6 /m- This implies that any known perturbation SU 
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can be included in when -jj- does not exceed the value given in table 8 
below : 





Earth 


Jupiter 


6U < e£0 2 
U ~ u 


2 1CT 8 


io~ 8 



Table 8. 



The quadrupolar contribution is bigger than 10 -3 , it cannot be consid- 
ered in known term, however higher multipole can be included in VE^ if 
the accuracy fi is smaller than 10 -6 for the Earth instead of 10 -4 and 
1CP 9 for Jupiter. The accuracy \i ~ 10~ 6 remains a very difficult challenge 
as far as the geometry of the set-up is concerned (A. Landragin, private 
communication) . 

With the same order of magnitude for fx, we obtain 



(ft)' 



= n 



t-LT 



'A - A A 



dt 



(k) 



(20) 



where VLlt is deduced from (THfl) and — — 

dt 



a -VU 



Mr, 



The satellites, such as the Moon for the Earth, can bring a contribution 
to U,(k)(j) X^X^ at the required level of accuracy but with such a value 
of (j., it could be included in *$>}.■ 

Of course when modelling one must be sure that every quantity is 
known at the required accuracy. This must hold for any value of /i, this 
is a necessary condition. Therfore the following relation must hold true : 





Earth 


Jupiter 




2 10- 12 


HT 12 



Table 9. 



Such an accuracy is not achieved for Jupiter. For the Earth, considering 
the difference between the various models (Godard Earth model 9 and 
10) it appears that the values of the high order multipoles are neither 
known nor consistent at the required level (10~ 12 ). One can hope that the 
lack of precision on the Jkn coefficients |l3| is not important for n / 
because the diurnal rotation modulates the frequency of the corresponding 
contribution in U,(k)(j) X^> X^\ However it is necessary to increase our 
knowledge of the axisymmetrical potential of the central body in order that 
£ fe AJfco < e£02 where AJ(. is the uncertainty on J^o — Jk- Such a relation 
holds true for k = 2. It could be presently achieved with low values of £ (on 
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high orbits) but the Lense-Thirring effect is proportional to £ 3 (see table 
4 above) and it seems impossible to measure the Lense-Thirring effect for 
£ << 1 in not to far a future. 

This question is crucial but Hyper might, itself, bring an answer to this 
question by the means of the time analysis. Now we forget these problems 
because in the simple case that we consider the quadratic quantities do not 
bring any contribution to the signal. 



We consider that the motion of the satellite takes place in the (x, ?/)-plane 
while the vector t lies in the (x, z)-plane. We assume that the eccentricity, 



A. Significant terms in ^ 



e, does not exceed £ 1 ^ 2 Oi. 
We define 



J = J x e j + J y e. 2 + J z e§ , n = cos 5j + sin 9 (21) 
— I = cos a ej + sin a eg , w = w x ej + w y eg + w z eg (22) 

J, I and w are constant vectors. The angle and the distance r depend on 



the time T. 



* 




FIG. 4: The satellite and the fixed star 



First we define the spatial triad u'. 



(n) ■ 
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Let us outline that we have defined f^(i)) — — (j^j + ^^ 2 0\. Then, 
in order to obtain the final tetrad u" a y we perform an arbitrary rotation 
around A : 

u(i) = -a* = *t' (1) + e 1/2 Oi 

U( 2 ) = u', 2 s cos a + u'ps sin cr + £ 1/,2 0i (24) 
U( 3 ) = — iZj 2 , sin a + it'^ cos cr + ^' 2 Oi 



where = ~ m is the angular velocity of the triad relatively to 



We can now assume that the experimental set-up is comoving with the 
triad ur n \ whose vector un) points towards the fixed star. 

During the flight of the atom, the quantity t/,(fe)(j) m equation 1191 does 
not remain constant because the position of the satellite changes. One can 
consider that the coordinate, X = of the atom is a function of the 

time: X = v g (T — To) . 

X 

Therefore we expand U,(k)(j) = U ,( k )(j) (To) + U,(k)(j)(t) ( T °) x _ 

v g 

where v^lim is the orbital velocity. 

Before performing explicit calculation we notice that 
U,(k)(j) (To) X^ k 'X^ J ' will not bring any contribution to the phase 
difference |3] because Os and 0' s are two centers of symmetry. For 

v (t) 

fi < 1(T 6 the term U, {k)ljm (T ) x X X^X^ ~ s 3 ^ 5/2 5 /r] can 

V 9 

be included in for Jupiter ( 5—5 < 1), but it must considered 

for the Earth I ' — ~ 5 > 1 I . However the quadrupole does not 



--Co,, 

bring any contribution to the phase difference that we calculate from ty u . 
In equation 1191 the spin zuu^ does not bring any contribution in the 

term v g because un\, X and v g are in the same plane. 
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Then, one obtains 




In expression (eq. Q25[0 of ^> u one can assume that r = r is a constant 
because we assume that the excentricity is small e < 0\ 1 therefore the 
corrections are included in 



Moreover, for the same reason, one can drop the terms of order 0\ in 
the expression of the tetrads. Therefore, it is clear that we can consider 
the space as the ordinary space of Newtonian physics and that the usual 
formulae to change the basis d % into e ^ or u r^are valid. 

In 1)25(1 . the terms of lines A and B are due to various rotations : respec- 
tively the Lense-Thirring rotation and the aberration. The term of line C 
is due to the displacement of the satellite during the flight time of the atom 
and the term of line D corresponds to some residual acceleration due to 
the fact that point O is not exactly in free fall. 



B. The phase differences 

We use the expression (12511 of ^> u in order to calculate 6<p given by [31 We 
find 
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where S = - — w„T,? is the area of the Sagnac loop. As we mentionncd 
Am 

before, J and -M@ are expressed in geometrical units. 

The two interferometers of the same ASU are assumed to lie in the same 
plane but not necessarily with their center of symmetry O5 and 0' s at 
the same point. Therefore adding and subtracting the phase differences 
delivered by the two interferometers one finds the two basic quantities which 

are measured by the set-up i.e. : /ii = — (Sip' — Sep) and ^2 — ^ l^v' + <V) • 
We define the shift 8 = X ' s - los arL d the acceleration ~a — 

— (floe + ao\) where ao s and Sq' are the accelerations at point Os and 

0' s . We drop several terms which can be included into ^ fe . Then we obtain 
the quantities which can be measured : 

Mi + ^M2 = ^(cT D f\^u {3 yn + n LT -u {2) \^ (27) 



c 



X 



2ir [cT D f M @ I _ 8 



Ar r V r 

2 ^ / nn ^2 (v g T D \ 2 M, 



W( 3 ) 3(n-U( 3 )j n 



A 




1 , x , , x , 4^(cT D ) 2 



M2 = - (V + = • « (3 ) } (28) 



C. Discussion 

We define a as the direction of the fixed star (fig. and the projection, 
J n of J on the plane of the orbit : 

J,, = J M (cos 6* j ifj + sin 9 j and iy M = w n (cos8 w ~e"j + sin^ ~~e 2 ) ■ 
Then 

111 + 2^ 2 = \- £L - x {Ko + K a + K 2a + K20 29 

c Ar 

+Kg_ a + Kg +a + K 2 e-a + K 2 g +a + K 2 g_ 2a + K 2 g +2a } 



[(1 — sin a) cos (a + 6j) — (1 + sin a) cos (cr — 9j)] x 



-4 cos a cos cr x — , 



3M @ 
4r 

_3M@ 
4r 



1 — sin 2 a) 



3M© . <5« 
cos a sin a sin ct x 

,5(2) 5 (3)' 

sin (2ct) x h cos (2cr) x 

r r Q 



(l - sin 2 a) cos (2(9) 



#3) 
^0 



= x {-2(1 + sin a) sin (9- a) 

r C 

H — — cos a sm{v — a — v w ) x — 

a i /-, • n • \ w 3 I 

+— (1 + sin a) sin (0 - cr) x — > 
Z c 



M ®v g 
r c 



x {-2 (1 - sin a) sin (9 + a) 



ai . . . w„ 

+— cos a sin [9 + a — 9 W ) x — 
2 c 

— — (1 — sin a) sin (9 + cr) x — 

£ C 



K 2 e-a 



-3^ (1 + sin a) cos (29 - a - 9j) x 
c r 

+ 3M© s in a) sin (26* - cr) x - — 



'o 



i^2e+ CT = 3-^ (1- sin a) cos (20 + a - 6j) x 



+ 3M© _ s in a) sin (29 + a) x - — 

4r r 



(M @ \ 3/2 cv g T3 
\ ro J 2r 2 



(1 — sin a) cos (20 + cr) 
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3M© 



(1 + sin a) 2 



(38) 



Sr 



K 2 6+2a = 



x <^ sin (29 + 2a)x h cos (29 + 2a) x 



x <^ sin (29 - 2a) x cos (29 - 2a) x 




5 (3) 



6 (3) 



} 



} 



(39) 



Each of these terms, except Kq, has a specific frequency. They can be 
measured and distinguished from each other. 

The Lense-Thirring effect due to the angular momentum of the central 
body appears in the terms K„ and K26± a while the possible existence of 
a preferred frame appears in Kg± a which depends on the components of 
a{w . 

The signal due to the Lense-Thirring effect is associated with the signal 
due to Today, it seems impossible to reduce 5^ significantly, this is 
the reason why it should be calculated from the Fourier analysis of the 
signal itself altogether with the velocity a{w . 

The interest of the spin is obvious. If a is constant (no spin) the signal is 
the sum of two periodic signals with frequency vo and 2vo where vo is the 
orbital frequency of the satellite ; therefore one ASU gives two informations 
(two functions of the time). When the satellite spins, we get 9 functions of 
the time t. The information is much more important in this case. 



In Hyper, the Lense-Thirring effect is associated with many perturba- 
tions which cannot be cancelled. We have exhibited the various terms 
that one needs to calculate in order to obtain the full signal and we have 
emphasized the necessity to increase our knowledge of the Newtonian grav- 
itational potential. This is still more crucial for Jupiter despite the fact 
that the Lense-Thirring effect is much bigger. 

Using four parameters, £, e, r\ and \i defined in table 2, we have also 
sketched a method to take into account the residual gravitational field in 
a nearly free falling satellite, namely the tidal and higher order effects. 

Compared with GPB, the principle of the measure is not the same, the 
difficulties are quite different but the job is not easier. For instance, consid- 
ering the quantities K a or K<ie±o above, one can check that, for an Earth 
satellite, 5^ must remain smaller than 2nm for the corresponding signal 
to remain smaller than the Lense-Thirring one. It does not seem that such 
a precision can be controlled in the construction of the experimental device 
itself. It is therefore necessary to measure with such an accuracy. 



V. 



CONCLUSION 
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What can be deduced from the time analysis depends on the accuracy 
of the various parameters. From K 2 g±2a we deduce a. Then from Kg± a 
we obtain aiw„/c and «iw 3 /c as two functions of W . Therefore one can 
check if a\w = or not. 

From K 2 0± a one can calculate 2 ® C ( 5( 1 ) as a function G(+) of J»/Tq, 8j 



f Mq\ 3/2 c 2 T 2 

and ^ ^ and as a different function, G(_) of the same argu- 

\ r o ) r t> _^ 
ments. One could check the equality Gi + \ = Gi_y If we assume that J@ 
is known, then 9j is known and from the equality G/ + \ = Gr_\ we deduce 

l + 7 + ai/4 / — > \ 
the value of J n . Using the relation J M = (^©J one could 

check whether j + a.i/4 = 1. 

K a would give as a function of J,,/?' 2 , 8j and jV r o- Using the 



0^9 



r o v a 



previous results, we obtain J 3 . The relation J 3 = (^®J 

gives an other test of the value of 7 + ai/4. 

But over all, the best test would be that the signal (as a function of the 
time) fits the theoretical prediction. 

As a final conclusion let us put forwards that the geometric scheme which 
has been used is just a preliminary contribution to the discussion on the 
feasibility of Hyper. Only a more powerful model can answer the question. 
This model should take into account all the gravitational perturbations 
that we have outlined here and it should consider the interaction between 
laser fields and matter waves in more a realistic manner. Such an approach 
has been recently developed 0, @ it could give definitive results in the 
future. 
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